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SELECTION  OF  OPENING  MODEL  FOR 
PARACHUTE  SCALING  STUDIES 


Introduction 

The  purpose  of  this  ongoing  research  project  has  been  to  identify  those  scaling  laws 
important  in  modeling  parachute  opening  dynamics,  and  to  determine  scaling  parameters 
that  will  provide  relationships  in  canopy  stiffness. 

Reference  (1)  summarizes  the  work  done  on  this  project  during  the  time  period  from 
December  15,  1988,  to  December  14,  1989.  During  this  period,  an  improved  scaling 
parameter  was  developed  for  determining  stiffness  effects  on  the  opening  behavior  of  flat 
circular  parachute  canopies.  This  scaling  parameter,  called  the  relative  stiffness 
index,  worked  well  for  correlating  opening  time  of  flat  circular  parachutes,  and  gave 
fair  correlation  for  predicting  opening  shock  for  these  canopies;  but  more  work  needed  to 
be  done  on  the  opening  shock  method.  In  addition,  the  relative  stiffness  index  concept 
should  also  be  extended  to  other  parachute  geometries.  The  concept  is  well  enough 
advanced,  however,  to  justify  its  incorporation  in  some  manner  into  an  opening  dynamics 
theory. 

In  the  present  work,  a  theory  was  chosen  and  computer  programmed  as  a  basis  for 
allowing  introduction  of  the  stiffness  concept  and  to  provide  a  basis  for  later 
correlation  of  the  test  data  available  at  Natick  to  further  validate  the  use  of  the 
relative  stiffness  index. 


Brief  retieM  of  araiiabie  theories 

A  listing  of  some  of  the  various  opening  dynamics  theories  can  be  found  in  [2-13].  A 
critical  review  of  some  of  the  opening  dynamics  theories  can  be  found  in  [14].  The  two 
theories  that  proved  most  useful  for  this  study  were  those  developed  by  Fu  [12]  and 
Lingard  (I3J. 

Fu  developed  a  theory  to  predict  parachute  opening  shock  and  time,  using  a  minimum  of 
experimental  inputs  based  on  as  little  test  data  as  possible.  Fu  considered  the  payload 
and  canopy  as  two  point  masses  connected  by  a  spring  with  damping.  The  canopy  shape  was 
represented  as  a  truncated  cone  topped  wirh  a  hemispherical  cap.  The  apparent  mass  of 
the  parachute  was  assumed  to  be  proportional  to  the  canopy  volume.  The  canopy  vent  size 
and  effective  porosity  of  the  canopy  fabric  were  considered  in  the  conservation  of  mass 
equations  used  to  model  canopy  filling.  Using  Newton's  Law',  Conservation  of  Mass,  and 
the  appropriate  geometric  relations,  the  writer  derived  a  system  of  nonlinear 
differential  equations.  Some  of  the  results  of  this  theory  are  presented  in  [1]. 

Lingard's  analysis  for  predicting  opening  shock  was  considerably  different  from  that 
of  Fu  but  led  to  similar  results.  His  method  was  intended  to  provide  a  relatively  simple 
theory  to  let  a  designer  perform  trade-offs  upon  the  effects  on  the  peak  load  of 
parameters  such  as  snatch  velocity,  suspended  mass,  parachute  size,  and  altitude 


and  angle  of  deployment.  The  main  assumptions  made  by  Lingard  were  as  follows: 

a.  The  parachute  inflates  in  a  constant  nondimensional  opening  time 
independent  of  snatch  velocity  and  mass  ratio. 

b.  For  a  given  parachute  design,  there  exists  a  unique  force  coefficient  vs. 
nondimensional  time  function  which  is  independent  of  snatch  velocity  or 
mass  ratio. 

These  two  assumptions  were  then  used  to  extract  empirical  data  from  a  small  number  of 
full-scale  tests  of  a  particular  solid  cloth  parachute  system.  With  these 
assumptions,  Newton’s  Second  Law  was  then  used  to  develop  a  model  that  allowed 
calculation  of  force  and  velocity  profiles  during  deployment.  The  effects  of  the 
above-mentioned  variables  on  deployment  could  then  be  studied. 

A  comparison  of  results  from  Lingard’s  and  Fu’s  methods  is  given  in  (1]  and  will 
not  be  repeated  here.  The  comparison  shows  similar  trends  in  results  from  the  two 
totally  different  approaches  and  lends  confidence  to  the  results  obtained.  Because 
of  assumption  (a)  used  above,  Lingard’s  method  is  not  useful  in  predicting  Froude 
number  and  mass  ratio  effects  on  opening  /mtc,  but  the  time  assumption  itself 
agrees  fairly  well  with  the  theoretical  results  obtained  by  Fu. 

More  complex  opening  dynamics  theories  are  also  available,  some  of  them  being 
briefly  summarized  in  (15].  However,  most  of  these  are  far  too  complex  for  general 
use  and  not  amenable  to  a  reasonable  simulation  of  opening  dynamics,  especially  when 
parametric  studies  need  to  be  conducted  for  a  wide  range  of  variables.  The  following 
quote,  taken  from  (16)  aptly  depicts  the  situation  regarding  many  of  these  prediction 
methods: 

The  difficulty  in  predicting  the  performance  of  a  parachute  lies  in  the  complex 
interaction  between  the  porous  canopy  and  the  surrounding  flow  field.  The 
parachute  inflation  process  involves  the  unsteady,  viscous,  compressible  flow 
about  a  porous  body  that  undergoes  large  shape  changes.  Moreover,  this  body  is 
composed  of  nonlinear  materials  with  complex  strain,  strain  rate,  and  hysteresis 
properties.  Thus  it  is  not  surprising  that  a  rigorous  analysis  of  the 
Navier-Stokes  equations  for  the  unsteady  flow  about  an  inflating  parachute 
presents  a  formidable  challenge  to  existing  computational  capabilities.  A  recent 
study  at  Sandia  National  Laboratories  concluded  that  an  axisymmetric  flow  field 
solution  for  a  typical  weapons  parachute  would  require  3300  hours  on  a  CDC  7600 
computer  or  330  hours  on  a  CRAY  system.  Hence,  there  exists  an  urgent  need  for  a 
dependable  intermediate  theory  useful  for  parachute  design. 

Continued  work  on  canopy  stiffness  scaling  laws  will  be  hampered  if  a  reasonable 
theory  is  not  available  that  can  be  modified  in  some  fashion  to  simulate  canopy 
stiffness.  Hence,  the  theories  were  examined  to  select  one  most  suitable  for  this 
task. 

The  theories  of  Lingard  and  Fu  referenced  earlier  seem  to  provide  a  reasonable 
compromise  in  accuracy  without  being  excessively  unwieldy.  However,  proper 
verification  of  these  parachute  opening  theories  will  eventually  require  a  wide  range 
of  tests.  The  scatter  in  test  data  creates  thi'  need.  Most  of  these  tests  would 
have  to  be  conducted  outside  the  practical  range  of  mass  ratios  and  Froude  numbers 
normally  used  in  real  parachute  applications. 


Conversely,  at  Froude  numbers  on  the  order  of  four  to  six  where  opening  shock  loads 
are  reasonable,  extrapolation  of  the  theories  of  Fu  and  Lingard  show  almost  no  effect  of 
typical  mass  ratio  changes  on  nondimensional  opening  force.  For  this  reason  further 
tests  need  to  be  conducted  at  high  Froude  numbers  over  a  wide  range  of  mass  ratios,  and 
the  theories  of  Fu  or  Lingard  must  be  used  to  calculate  expected  parachute  behavior  at 
lower  Froude  numbers  than  they  present  in  their  works. 


Basis  for  final  selection  of  theory 

The  theory  selected  for  programming  must  ultimately  be  one  that  appears  suitable  for 
modification  to  simulate  canopy  stiffness  effects,  be  relatively  complete  (not 
oversimplified),  yet  amenable  to  computer  programming  and  calculations  at  reasonable 
cost.  The  theory  developed  by  Fu  appeared  to  be  the  best  for  this  purpose,  if  canopy 
stiffness  effects  can  eventually  be  simulated  using  various  input  values  for  initial 
semi-aperture  angle,  bo  (see  Fig.  1),  or  some  similar  means.  Heinrich  and  Hektner’s 
stiffness  index  [17]  used  the  inverted  hang  test  of  the  canopy  and  an  experimental 
measurement  of  Dmax  to  represent  canopy  stiffness.  The  Dmax  and  stiffness  index  values 
should  be  reflected  in  the  initial  semi-apertural  angle  bo  in  some  manner.  Fu’s  theory 
also  allows  part  of  the  areal  density  effect  of  the  fabric  to  be  reflected  by  making  an 
input  change  to  canopy  mass.  Some  of  the  effect  of  canopy  stiffness  could  be  accounted 
for  by  an  adjustment  to  suspension  line  elasticity  input,  but  this  does  not  duplicate 
the  main  effect  of  canopy  stiffness. 


Major  assumptions  in  the  theory 

The  complete  details  and  assumptions  used  by  Fu  will  not  be  given  in  this  report. 

The  reader  is  referred  to  [12)  for  this  information.  The  major  assumptions  and  the 
models  used  to  simulate  the  parachute  opening  geometry  will  be  summarized,  however.  The 
parachute  opening  is  broken  down  into  a  "first  filling  phase,"  "second  filling  phase," 
and  "transition  phase."  The  assumptions  are  as  follows; 

a.  First  filling  phase.  During  the  first  filling  phase,  the  parachute 

geometry  is  assumed  as  shown  in  Fig.  1. 

1.  The  initial  semi-apertural  angle,  bo,  is  assumed  to  be  a  constant  during 
this  phase,  and  is  obtained  empirically. 

2.  The  canopy  shape  starts  off  as  a  small  hemisphere  and  fills  until  the 
distance,  X,  is  equal  to  the  constructed  diameter  of  the  parachute.  This 
determines  the  end  of  the  first  filling  phase. 

3.  During  the  first  filling  phase,  no  porosity  or  vent  flow  is  assumed. 

4.  The  inflow  area  at  the  mouth  of  the  parachute  is  taken  to  be  a  fraction  of 
the  actual  area,  using  an  empirical  "contraction  factor,"  kj.  The 
canopy  velocity  is  used  for  the  inflow  velocity. 

5.  A  constant  parachute  drag  coefficient  (typically  0.5)  is  assumed  during 
phase  I ,  and  the  drag  changes  with  changes  in  projected  frontal  area  and 
canopy  velocity. 

6.  Each  suspension  line  is  represented  by  a  spring  and  damper,  allowing  for 
line  stretch  and  a  relative  velocity  between  canopy  and  payload. 
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7.  The  initial  suspension  line  stretch  is  determined  from  an  input  value  for 
canopy  snatch  force. 

8.  In  this  phase,  and  all  phases,  the  air  is  assumed  incompressible. 

As  mentioned  earlier,  once  the  "wrap-around  distance,"  X,  equals  the  canopy 
constructed  diameter,  the  first  filling  phase  ends  and  the  second  filling  phase 
begins. 

b.  Second  filling  phase.  Fig.  2  shows  the  model  used  to  represent  the  canopy 
during  the  second  filling  phase.  The  assumptions  used  to  represent  this 
portion  of  filling  are  as  follows: 

1.  The  shape  of  the  canopy  is  assumed  to  be  a  truncated  cone  with  a 
hemispherical  cap.  The  wrap-around  distance  remains  constant  and  equal  to 
the  constructed  diameter,  DC. 

2.  As  air  flows  into  the  parachute  mouth,  the  suspension  line  semi-apertural 
angle,  /?,  increases,  causing  the  hemispherical  cap  to  grow  larger  in 
diameter  while  the  sides  of  the  truncated  cone  portion,  h,  decrease  in 
length. 

3.  Conservation  of  mass  applied  to  the  inflow  air  governs  the  rate  of  canopy 
inflation.  Air  flow  is  allowed  out  through  the  vent  and  through  the 
porous  hemispherical  portion  of  the  canopy.  An  empirical  contraction 
factor,  k2,  different  from  k|  in  a.4.,  is  used  for  this  portion  of  the 
filling. 

4.  A  new  parachute  drag  coefficient  (typically  1.3)  is  assumed  during  phase 
2.  Again,  canopy  drag  changes  with  projected  frontal  area  and  canopy 
velocity. 

5.  Suspension  lines  are  again  represented  by  springs  and  dampers  allowing  for 
line  stretch  and  a  relative  velocity  between  canopy  and  payload. 

6.  The  second  filling  phase  ends  when  the  length  of  the  sides  of  the 
truncated  cone,  h,  has  decreased  to  a  value  of  0.03  D^.  This  is  an 
empirical  value  obtained  from  motion  pictures  of  parachute  opening 
analyzed  by  Fu. 

c.  Transition  phase.  The  transition  phase  is  that  phase  that  occurs  between 
the  end  of  the  second  filling  phase  and  the  time  that  steady  state  descent  is 
reached.  Fig.  3  shows  the  idealized  geometry  for  the  transition 

phase.  The  assumptions  used  in  this  phase  are  as  follows: 

1.  The  shape  and  size  of  the  canopy  remain  constant  at  the  steady  state 
values  hjp  D^,^2’  ’'sf  ^st  shown  in  Fig.  3.  These 

values  are  the  corresponding  values  that  existed  at  the  end  of  the  second 
filling  phase  as  defined  by  h^j  =  h  =  0.03  D^. 

2.  The  suspension  lines  are  still  represented  by  springs  and  dampers  allowing 
for  line  stretch  and  relative  velocity  between  canopy  and  payload.  Thus, 
although  the  canopy  shape  remains  unchanged,  the  semi-apertural  angle, 

varies  through  a  damped  oscillation. 
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Fig.  3.  Geometric  Model  of  the  Canopy  for  the  Transition  Phase  [12], 
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3.  Air  flow  is  again  accounted  for  through  the  vent  and  porous  hemispherical 
portion  of  the  canopy. 

4.  The  drag  coefficient  is  held  constant  at  the  value  assumed  for  the  second 
filling  phase. 

The  transition  phase  ends  when  all  geometric  and  trajectory  parameters  reach 
steady  state  values,  i.e.,  the  suspension  line  length,  1,  and  semi-apertural 
angle,  /S,  become  constant,  and  the  terminal  velocity  and  vertical  descent 
path  are  reached.  Theoretically,  this  would  occur  asymptotically  over  a  long 
period  of  time,  but  practically  speaking,  steady  state  values  are  approached 
closely  in  a  period  of  several  seconds.  In  any  case,  maximum  filling  force  and 
opening  time  have  already  been  calculated  near  the  end  of  the  second  filling 
phase  in  most  situations. 

Fig.  4  is  a  simplified  illustration  of  the  velocity  and  acceleration  components  for  the 
combined  canopy  and  payload  system.  Fig.  5  shows  the  forces  acting  on  the 
parachute-payload  system. 


Fig.  4.  Velocity  and  Acceleration  Components  of  the  Elastically 

Connected  Canopy  (m^^-)  and  Payload  (m^)  Mass  Points  in  the 
Body  Coordinate  System  (12). 
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Computer  program  development 


Each  phase  of  the  opening  process  involves  the  numerical  solution  of  between  five 
and  eight  simultaneous  nonlinear  differential  equations,  together  with  numerous 
geometric  constraint  equations. 

To  properly  program  the  equations  developed  in  [12],  it  is  necessary  to  become 
thoroughly  familiar  with  all  the  equation  derivations;  verifying  assumptions,  checking 
the  algebra  of  the  derivations,  checking  for  typographical  errors  in  the  preparation  of 
[12],  and  then  coding  the  program  into  a  computer  language  (in  this  case,  Fortran  was 
used).  Because  of  the  complexity  of  this  task,  and  due  to  time  limitations,  it  was  not 
possible  to  do  these  steps  for  all  three  phases,  and  some  compromises  had  to  be  made. 
These  were  as  follows. 

The  first  filling  phase  equations  were  checked  carefully,  some  typographical  errors 
were  noted  and  corrected,  and  the  equations  were  programmed.  Because  of  time 
constraints,  only  some  of  the  equations  were  checked  for  the  second  filling  phase  and 
transition  phase.  The  major  part  of  the  effort  went  into  programming  these  equations 
and  running  them  to  check  if  the  results  seemed  reasonable. 

One  point  should  be  made  here.  As  the  program  goes  from  one  filling  phase  to 
another,  the  relative  velocity  between  the  canopy  and  payload  must  be  recalculated 
between  phases  and  a  small  discontinuity  in  this  value  introduced.  This  discontinuity 
is  necessary  to  satisfy  force  and  length  continuity  in  the  suspension  lines,  between 
phases,  due  to  the  discontinuity  that  inevitably  must  occur  when  the  geometric  model  is 
abruptly  changed  between  phases.  In  most  cases,  the  change  in  relative  velocity  is  less 
than  a  few  percent  and  should  be  of  no  consequence  within  the  accuracy  of  the  theory. 
This  approximation  is  inherent  to  Fu’s  theory. 


Current  Liniiuilions  in  the  Pr()<;ram 

There  are  two  limitations  in  the  computer  program  that  should  be  removed  at  a  later 
date.  For  certain  calculated  drop  conditions,  the  calculated  force  in  the  suspension 
lines  goes  to  a  small  negative  value  for  a  short  period  of  time  during  the  transition 
phase.  During  the  transition  phase,  the  parachute-payload  is  behaving  primarily  as  a 
spring-mass-damper  system.  Conceivably  this  behavior  could  occur  during  the  other 
phases  for  some  drop  conditions.  Fu  made  no  mention  of  how  this  should  be  handled  in 
[12].  Obviously,  this  must  be  modified  by  program  statements  that  set  the  suspension 
line  force  equal  to  zero  if  a  negative  force  is  calculated,  since  the  suspension  lines 
cannot  physically  tolerate  a  compressive  force.  In  [5],  Payne’s  theory  also  considers 
suspension  line  stiffness,  but  not  canopy  stiffness,  and  the  suspension  line  force  also 
goes  to  zero  several  times  during  some  drops  due  to  the  spring-mass-damper 
characteristics  of  the  system. 

The  other  change  that  needs  to  be  made  to  the  program  is  an  improvement  in  the  accuracy 
of  the  integration  routine  used.  The  current  program  uses  Euler’s  method  for  solving 
the  differential  equations.  This  is  adequate  as  a  first  approximation  to  the  theory, 
but  needs  to  be  upgraded  before  the  program  is  used  extensively.  A  higher  order 
integration  routine,  such  as  fourth-order  Runge-Kuita,  should  be  incorpoialcu. 
Additionally,  time  should  be  spent  checking  the  accuracy  of  the  remaining  equations  in 
[12]  to  ensure  that  no  additional  typographical  or  derivation  errors  exist. 
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Program  Listing  and  Typical  Output 


Appendix  A  contains  the  Fortran  listing  for  the  program,  temporarily  called  FUDROP 
here,  with  input  for  a  typical  case.  The  input  is  for  a  C-9  type  parachute  canopy,  and 
the  parachute  geometry  and  drop  conditions  are  given  in  Table  1.  The  values  in  Table  I 
correspond  to  Example  1  from  [12]  and  are  used  to  compare  this  program’s  output  with 
Fu’s  results. 

Appendix  B  contains  the  program  output  for  this  case.  The  first  part  of  the  program 
is  an  echo  of  the  input  for  parachute  geometry  and  snatch  conditions.  This  is  followed 
by  the  output  for  the  first  filling  phase,  lasting  for  0.20  seconds,  followed  by  a 
statement  indicating  the  end  of  the  first  filling  phase. 


Table  1.  Parachute  Geometry  and  Drop  Conditions,  FUDROP 
(See  Appendix  A) 


Variable 

Symbol 

Value 

Parachute  nominal  diameter 

Do 

28.22  ft 

Constructed  diameter 

Dc 

28.45  ft 

Payload  radius 

R 

0.755  ft 

Vent  radius 

Rv 

0.820  ft 

Suspension  line  length 

Is 

28.22  ft 

Air  density  (1000  ft  alt) 

P 

0.00231  sl/ft3 

Number  of  suspension  lines 

n 

28 

Drag  coefficient  of  payload 

Cdl 

0.64 

Snatch  velocity 

Vs 

185  fps 

Total  sveight 

Wt 

296  lbs 

Payload  weight 

Wl 

280  lbs 

Initial  trajectory  angle  (below  horizontal) 

Bo 

50“ 

Snatch  force 

Fs 

662  lbs 

Effective  porosity  of  fabric 

Ccff 

0.1 

The  output  format  for  the  second  filling  phase  is  presented  in  a  slightly  different 
arrangement,  and  ends  at  0.90  seconds  for  this  example.  The  output  for  the  transition 
phase  is  again  in  a  different  format,  and  calculations  for  this  case  were  stopped  at  5.9 
seconds,  as  the  terminal  velocity  had  essentially  been  reached  and  the  trajectory  was 
approximately  vertical  (6  =  85’). 


Comparison  with  Fu's  Results 


The  calculated  results  from  Appendix  B  are  nondimensionalized  and  compared  with  Fu’s 
results  in  Fig.  6.  In  this  figure,  the  solid  lines  are  Fu’s  calculated  results,  the 
short  dashed  lines  drawn  through  squares,  triangles  or  diamond-shaped  symbols  are  the 
Appendix  B  results,  and  the  long  dashed  lines  are  experimental  data  from  [18],  The 
letters  along  the  abscissa  of  Fig.  6  correspond  to  comparisons  of  the  calculated  and 
measured  parachute  shapes,  semi-apertural  angles,  and  nondimensional  times  in  Fig.  7. 

Figure  7  shows  fair  to  good  agreement  between  theory  and  experiment.  The  major 
difference  is  that  the  parachute  changes  shape  somewhat  at  the  beginning  of  the 
transition  phase,  including  the  well  known  "overexpansion."  The  theory  does  not  account 
for  this,  nor  for  the  change  in  suspension  line  angle  during  the  first  filling  phase. 

Comparing  the  results  in  Fig.  6  for  nondimensional  radius,  T,  the  Appendix  B  results 
differ  from  Fu’s  somewhat,  but  agree  better  with  experiment.  The  present  theory  shows 
the  parachute  growing  in  size  more  slowly  than  Fu  predicts.  Comparing  nondimensional 
velocity  (vj_)  histories,  Fu’s  theory  agrees  better  with  experiment  than  the  Appendix  B 
results.  The  slower  theoretical  parachute  growth  rate  of  the  present  theory  results  in 
less  drag  and  a  resulting  slower  velocity  decay  than  predicted  by  Fu.  For  the  force 
history  comparison,  Fu’s  results  agree  better  with  experiment  during  the  early  part  of 
the  opening  time  period,  but  the  Appendix  B  results  agree  better  with  experiment  near 
the  peak  period  of  the  force  history.  Generally,  the  trends  shown  for  both  theories  and 
the  experiment  agree  well. 

It  is  expected  that  better  agreement  between  the  present  programming  of  Fu’s  theory 
(Appendix  A)  and  Fu’s  results  would  be  obtained  if  a  more  accurate  integration  routine 
were  used  in  Appendix  A.  The  type  of  deviation  between  the  two  theories  shown  in  Fig.  6 
is  typical  of  the  difference  obtained  when  using  Euler’s  method  for  integration  as 
opposed  to  higher  order  methods. 


Suggestions  for  future  work 

The  present  theory  should  be  improved  by  updating  the  integration  routine  and 
incorporating  program  statements  to  set  riser  force  equal  to  zero  whenever  the  program 
calculates  a  negative  suspension  line  force  (this  has  occurred  for  the  present  theory 
only  for  short  periods  of  time  in  the  transition  phase  for  cases  examined  thus  far). 

The  remaining  equations  developed  by  Fu  should  be  checked  for  accuracy. 

Fu  applied  this  theory  only  to  aircraft  drop  test  data  of  full-scale  28  ft  C-9 
parachutes.  The  program  should  be  used  to  compare  predictions  with  Lee’s  drop  test  data 
[19,20]  for  the  l/2-scale  and  1/4-scale  C-9  models,  and  to  the  full-scale  C-9  parachutes 
that  were  drop  tested  with  the  heavier  fabrics  as  part  of  Natick’s  scaling  studies. 

Finally,  the  program  should  be  modified  in  some  fashion  to  simulate  canopy  stiffness 
effects  so  the  program  can  be  better  used  to  study  the  scaling  problems  that  arise  when 
using  model  wind  tunnel  or  drop  test  data  to  predict  prototype  parachute  performance. 
The  program  could  then  be  used  as  a  basis  for  further  theoretical  verification  of  the 
relative  stiffness  index  proposed  in  [21]. 
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Fig  7.  Calculated  and  Measured  Geometric  Shape  of  the  Canopy  at  Various 
Nondimensional  Times  (Measurements:  [18],  Figure;  (I2J). 


Conclusions 


Any  evaluation  of  experimental  data  on  parachute  scaling  stiffness  effects  is 
somewhat  hampered  if  a  simplified  theory  is  not  available  to  guide  in  interpreting  test 
data.  This  report  has  presented  a  theory  that  can  be  used  for  this  purpose  with 
suitable  modifications. 

The  theory  has  been  programmed  in  Fortran  and  checked  against  Fu’s  results  for  a 
single  case.  More  work  needs  to  be  done  to  check  the  equation  derivations  and  upgrade 
the  integration  routines. 

After  the  improvements  are  made,  the  theory  should  be  compared  with  available  Natick 
drop  test  data  on  scale  model  parachutes,  and  the  effects  of  canopy  stiffness  should  be 
incorporated  into  the  theory. 
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Appendix  A 

FORTRAN  LISTING  OF  FUDROP  PROGRAM 
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PROGRMI  rODROP 

C  THIS  PROGRMC  CALCUIATBS  THE  OPENXMG  DICKMCICS  OP  A  PARACHOTB 
C  USING  lU'S  THXOIUr 

C  THIS  PORTICHf  OP  THE  PROGRAM  HTCXiUDES  GEOMETRIC  PARACHUTE  INPUTS 
C  AMD  INITIAL  DEPLOYMENT  CONDITIONS 

CP«20.6 
CEP-O.l 
CDl-0.5 
RPaO.755 
CD2-1.315 
CDL»0 . 64 
DS«1.37 
DO>28.22 
0028 . 45 
PSb662 . 

G-32.19 

XKs2. 

XXl-0 . 49 

XLSa>28.22 

XNb28. 

RV^O . 82 
Ml^280. 

NT>296. 

VS«185 
ZLolOOO . 

BO>0.04 
P-3.1416 
THO-0 . 8727 
RN^O. 

RHO-0 . 00231 
XONO-0. 

PRINT*, 'THE  DEPAULT  INPUTS  POR  THIS  PROGRAM  ARE  AS  FOLLOWS' 
PRINT*, 'PARACHUTE  DIAMETER  IS',DO, 'PT' 

PRINT*, 'CONSTRUCTED  DIA  OP  PARACHUTE  IS',DC,'PT' 

PRINT*, 'NO.  OP  SUSPENSION  LIMES  IS',XN 
PRINT*, 'SUSPENSION  LIME  ZAIGTH  IS',XLS,'PT' 

PRINT*, 'APERTURE  RADIUS  IS',RV, 'PT' 

PRINT*, 'EPPECTIVS  POROSITY  IS', CEP 

PRINT*, 'SPRING  CONSTANT  OP  A  SUSPEMSI(»(  LINE  IS' ,CP, 'ZA/PT' 
PRINT*, 'DAMPING  COEP.  OP  A  SUSP.  LINE  IS' ,DS, 'LB-SEC/PT' 

PRINT*, 'PIRST  PHASE  PARA.  DRAG  COEPP  IS',CD1 
PRINT*, 'SECOND  PHASE  PARA.  DRAG  COKFF  IS' ,CD2 
PRINT*, 'TOTAL  CANOPY  S  PAYLOAD  HEIGHT  IS', NT, 'LBS' 

PRINT*, 'APP.  MASS  PROPORTIONALITY  INSTANT  IS',XK 
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PRINT*, 'riRST  PHASE  CONTRACTION  FACTOR  IS',XK1 
PRINT*, ' INITIAL  SEMI -APERTURE  ANGLE  IS' , BO, 'RADS' 
PRINT*, 'PAYLOAD  HEIGHT  IS' ,NL, 'LBS' 

PRINT*, 'PAYLOAD  RADIUS  IS',RP,'rT' 

PRINT*, 'PAYLOAD  DRAG  COETT  IS',CDL 
PRINT*, '  ' 

PRINT*, 'PARACHUTE  SNATCH  CONDITIONS  ARE  AS  FOLLOWS:' 
PRINT*, ' SNATCH  VELOCITY  IS' , VS, 'FT/SEC' 

PRINT*, ' INITIAL  FLIGHT  PATH  ANGLE  THETA  IS ' , THO, ' RADS ' 
PRINT*, 'X  COORD  OF  PAYLOAD  AT  SNATCH  IS', XL, 'FT' 

PRINT*, 'ALTITUDE  OF  PAYLOAD  AT  SNATCH  IS',ZL, 'FT' 
PRINT*, 'OPENING  SHOCK  OF  PAYLOAD  AT  SNATCH  IS',FS,'LBS' 
PRINT*, 'REL.  VEL.  BETWEEN  CANOPY  AND  LOAD  IS',RH, 'FPS' 
PRINT*, 'INITIAL  ANG.  ROTATION  RATE  IS' ,XC»«0, 'RAD/SEC' 

C  INITIAL  VALUES  OF  VARIABLES  FOR  FIRST  FILLING  PHASE 
DTe.OOl 
DTPR?:.01 
TPRxDTPR 
VLT^S 
VCN*0. 

THaTHO 

XOKaXOMO 

XLO« (FS/ (XN*CF*COS (BO) ) ) +XLS 

ROa ( (FS/ (XN*CF*COS (BO) ) ) +XLS) *TAN (BO) 

VO- ( (P/3 . ) *RO**3 . ) * (4 . - (2 . +SIN (BO) ) * (1 . -SIN (BO) ) **2 . ) 

RaRO 

H«0. 

VaVO 
XLaXLO 
XML>4n./G 
XMCa(NT-NL)  /G 
CDaCDl 
CBOaCOS (BO) 

SBOaSIN(BO) 

TBOaTAN(BO} 

CTBOal , /TBO 

500  DPR^  (P*R**2  . )  *  (4 . -  (2  .  -t-SBO)  *  (1 . -SBO)  **2 . )  +  (P*H*CBO**3  . )  * 

C  (2 . *R-H*TBO) 

DPH» (P*CBO**3 . ) * (R-H*TBO) **2 . 

VDa (XKl*P*CBO**2 . ) * (VLT-RW) * (R-H*TBO) **2 . 

RDa(  (RW/CBO)  *DPH-»-VD)  /  (CTBO*DPH-(-DPR) 
BDa(VD*CTBO-RW*DPR/CBO) / (CTBO*DPR+DPR) 

XU>aRD*CTBO-HD 

THDaX<B( 

DVaVLT-RW 
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RCHi4l*CTB0-B 

VLTD^VCN*XOM+ (X0M**2 . ) *RCH*CBO- (XH*CBO/XML) * (CT* (BCB 
C  -XLS)  4-DS*RM/CBO)  •ra*SXN  (TH)  -  ( .  5*CDt.*RBO*VLT/3afL)  *SQRT  ( 

C  VLT**2 .  +  (VCH-l-XOM*RCB*CBO)  **2 . )  *9*9»**2 . 

IWDiiVlTD-Vai*XOM+  (  (XK*FHO*OV*VD)  /  (XMC-fXK*RHO*V) )  -  (  (XM* 

C  CBO*  (Cr*  (RCH-XLS)  4-OS*Rir/CBO}  •fXMC*G*SIN  («B)  - .  5*CD*RBO 

C  *DV*SQI«P  (DV**2 .  +VCM**2 . )  *P*R**2 . )  /  (XMC-l-XK*RHO*V)  ) 

VGMD« 

C  -DV*XOM>  (XK*RHO*VCM*VD)  /  (XMC4XK*ItH0*V)  *  (XMC*6*C0S  ( 

C  TH) - . 5*CD*RHO*VCH*SQRT (DV**2 . +VCN**2 . ) *P*R**2 . ) / (XNC 

C  +XK*RHO*V) 

XOND-  (-VCND-  (VLT-fRH)  *XOM+G*COS  (TH)  -  (0 . 5*a>L*RBO/XML} 

C  *  (VCN-t-XOM*RCB*CBO)  *SQRT  (VLT**2 .  *  (VCM+XOM*RCH*CSO)  ** 

C  2 . ) *P*RP**2 . ) / (CBO*RCH) 

XZ.>R*CTBO-R 

XU>«HD*CTBO-BD 

r«XM*CBO*  (Cr*  (XL-XLS)  H-DS*XXi)) 

600  Xr(T-TPR)  610,630,630 
610  CONTIMDX 
620  GO  TO  660 
630  PRnrr*,'  ' 

PRINT*, 'TIMB  RATE  OF  CHANGE  OF  IS' , VD, 'CFS/S' 

PRINT*, 'TIME  RATE  OF  CHANGE  OF  RADIUS  IS' ,80, 'FT/S' 

PRINT*, 'TIME  RATE  OF  CBANQE  OF  B  IS ' , HD, ' FT/5 ' 

PRINT*, 'PAYLOAD  TANGENTIAL  ACC.  IS' ,V1.TD, 'FPS/S' 

PRINT*, 'ACCEL  OF  CANOPY  NRT  PLD  IS' ,8110, 'FPS/S' 

PRINT*, 'CANOPY  N08MAL  ACC.  IS' ,VCND, 'FPS/S' 

PRINT*, 'RISER  FORCE  ACTING  ON  PAYLOAD  IS', F, 'LBS' 

PRINT*,'  ' 

PRINT*,'  TIME  (SECS)  ','  VOL{COFT)  ',  '  RADIUS  (FT)', 

C  '  HEIGHT  (FT) ', '  LEl^TR(FT)' 

PRINT*, T,V,R,H,XL 

PRINT*,'  THETA (DEG) ', '  VLT(FPS)','  RW(FPS)', 

C  '  VCN  (FPS)  ','  OMEGA (RAD/SEC)' 

THs57.296*TH 

PRINT*, TH, VLT, RN, VCN, XOM 
PRINT*,  '  ' 

TH>TH/57.296 
650  TPR^TPR-I-DTPR 
660  TEST-H-f  (P/2 .  -t-BO)  *R 

IF  (TEST. GE. (DC/2.))  GO  TO  900 

VtaV+VD*OT 

IMl-»-RD*DT 

HaB-fBD*OT 

XL-XL4XLD*DT 

TH-TH-fTHD*DT 
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Vl,T»VI.T+VLro*DT 

RNtaRN-t-Rfa>*DT 

VC»^VCN+VCND*DT 

XOM>XOM+XCMD*DT 

THDsXOM 

T^+DT 

ir(T.GB.0.3)  GO  TO  1000 
GO  TO  500 

900  PRINT*, 'END  OF  FIRST  FILLING  PHASE' 

1000  CONTINUE 

C  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  PARACHUTE 
C  BEHAVIOR  DURING  THE  SECOND  FILLING  PHASE 
PRINT*,'  ' 

PRINT*, 'PARACHUTE  DYNAMICS  DURING  SECOND  FILLING  PHASE' 

DT=.001 
DTPRfO . 05 
TPR=T+DTPR 
CD=CD2 
BbBO 

C  PROGRAM  REENTERED  HERE  EACH  TIME  WITH  NEW  VALUE  OF  BETA 
1500  CBsCOS(B) 

SBsSIN(B} 

TBaTAN(B) 

CTBal . /TB 
GE«0 . 5*P+B 

CTO^O .  5*DC-  < .  5*P-»-B)  *R 
GEC»t=R*CTB-GBO 
A1>CB*  (CTB+<5E) 

A2» (P*R**2 . ) * (4 . - (2 .+SB) * (1 . -SB) **2 . ) + ( (P*CB**3 . ) /3 . ) * ( (-GE) *(0.75 

*  *R**2  .+  (1 . 5*R-TB*GEO)  **2  . )  +<3EO*  (1 . 5*R+2 . 

*  *(1.5*R-TB*GEO)*(1.5+TB*GE))) 

Bl— (R*CB*CTB**2 .+SB* (R*CTB-GEO) ) 

B2»P* (R*CB) **3 . + ( (2 . *P*CB**3 . ) /3 . ) *GEO* (1 . 5*R-TB*GEO) * (R*TB- 
+  (GEO/CB**2 . ) ) - ( {P*CB**3 . ) /3 . ) * (R+3 . *TB*GEO) * (0 . 75*R**2 . + 

+  (1.5*R-TB*GEO)**2.) 

HaGEO 

IF(B.GT.0.1*DO)  XK2al.O 
IF(B.LE.0.1*DO)  XK2=0 . 1+9 . *H/DO 
IF(H.LE. .03*DO)  GO  TO  2000 
IF(B.GT.BO)  GO  TO  1540 

C  THIS  PORTION  OF  PROGRAM  ITERATES  ON  VALUE  OF  RN 

C  NEEDED  FOR  TRANSITION  FROM  PHASE  I  TO  PHASE  II 

1530  PRINT*, 'ITERATION  SEQUENCE  FOR  RN  BETNEEN  PRASES  I  AND  II' 
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PRINT*, 'rZNJO.  PHASI  I  RN  BQUJILS '  , RN,  ' PPS 
COMTXNDB 


1540  VD-P*  (VLT-RM)  *XK2*  { ( (R-6E0*TB)  *CB)  **2 .  -RV**2 .  -2 .  *CaBF*R*R*  (SB+ 

+  SQRT (1 . - (RV/R) **2 . ) ) ) 

RD-  (RII*B2-B1*VD)  /  (A1*B2-B1*A2) 

BDm  (A1*VD-RN*A2)  /  (A1*B2-B1*A2) 

IP (BBS (B-BO) .GB. .00002)  GO  TO  15€0 
TBNPU)>RD*  (0 . 5*P-»-B4CTB)  -R*BD*CTB**2 . 
ir(ABS(XlJ3-TEMPll>)  .LE.  .001)  GO  TO  1550 
R1»bRN-BD*XL*SB 
GO  TO  1540 
1550  CONTINDB 

PRINT*, 'INITIAL  PHASE  II  RN  EQUALS' ,RN, 'TPS' 

XXJXTEMPLD 

C  END  or  ITERATION  ROUTINE  TOR  RN  BETNEEN  PHASES  I  AND  II 
1560  THD-XOM 

VLTD^VCN*X(»(-t-  (GEOM)  *CB*X0M**2 . -  (  (XN*CB)  /XKL)  *  (CT*  (GBOM-XLS) 

•I-  ■¥  <DS*RN/CB)  -fDS*  (GE(B1)  *BD*TB)  •K3*SIN  (TH)  -  ( (CDL*RHO*VLT)  /  (2  .  *XML) ) 

+  *  (SQRT  (VLT**2 .  +  (VCXt+XOH*GBOH*CB)  **2 . )  )  *P*RP**2 . 

DVWVLT-RN 

RND*VLTD-VCH*XOM+  (XK*RHO*DV*VD)  /  (XMC+XK*RHO*V)  -  (  {XN*CB)  *  (CP* 

*  (GEOM-XLS)  -«'DS*BN/CB-fDS*GE(»<*BD*TB)  •fXNC*6*SIN  (TH)  >0 . 5*CD*RBO 

+  *DV*  (SQRT  (DV**2  .  +VCN**2 . )  )  *P*R**2 . )  /  (XMC+XK*RHO*V) 

VCND—DV*XOM-  (XK*RHO*VCN*VD)  /  (XMC-»^XK*RHO*V)  +  (XMC*G*COS  (TH) 

+  -0 . 5*CD*RHO*VCN*  (SQRT  (DV**2 .  +VCN**2 . )  )  *P*R**2 . )  /  (XMC-^XK*RBO*V) 

XOMD- (-VCND- (VLT+RN) *XOM4G*COS (TH) - (0 . 5*CDL*RHO/XML) * (VCN+XOM 
+  *GBOI€*CB) * (SQRT (VLT**2 . + (VCN+X0M*GB0M*CB) **2 . ) ) *P*RP**2 . ) / ( 

■f  CB*GBON) 

r>XN*CB*  (CP*  (R*CTB-H-XLS)  +DS*  (RD*  (0 . 5*P-I-B4CTB)  >R*BD*CTB**2  . )  ) 

1600  ir(T-TPR)  1610,1630,1630 
1610  CONTINUX 
1620  GO  TO  1660 
1630  PRINT*,'  ' 

PRINT*,'  TIME (SECS)  ','  BETA (DEG)  ',  '  VOL (COPT)',' 

•f  RADIUS  (PT)',  '  rORCE(LBS)' 

B-57.296*B 
PRINT*,  T,B,V,R,r 
B-B/57.296 

PRINT*,'  THETA(DEG)','  VLT(PPS)','  RN(rPS) ' , 

*  '  VCN(rPS)  ','  OMEGA (RAD/SEC) ' 

TH>57.296*TH 

PRINT* , TH, VLT, RN, VCN, XON 
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PRIKT*, '  ' 

1650  TPI^^TPR+DTPR 
TH-TH/57.296 
1660  V^+VD*DT 
R«R+RD*DT 
B»B+BD*DT 
TH»TH+THD*DT 
VLT*VLT+V1,TD  *DT 
Rl»sRN+RWD*DT 

vcN=vcN+vcaro  *dt 

Xa«=XOM-l>X(»«D*DT 

TssT+DT 

IF(T.GE.5.)  GO  TO  2500 
GO  TO  1500 

2000  PRINT*, 'END  OF  SECOND  FILLING  PHASE' 

C  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  PARACHUTE  BEHAVIOR 
C  DURING  THE  TRANSITION  PHASE 
PRINT*,'  ' 

PRINT*, 'BEGINNING  OF  TRANSITION  PHASE' 

PRINT*, •  ' 

DT*.001 

DTPR?b0 . 1 

TPR*T+DTPR 

RST*R 

HST»H 

BSTsB 

VST«V 

PRINT*, 'STEAD!  STATE  CANOPY  VOLUME  EQUALS' , VST, 'CU  FT' 
PRINT*, 'STEADY  STATE  CANOPY  RADIUS  EQUALS' ,RST, 'FT' 

PRINT*, 'HST  EQUALS' ,HST, 'FT.  AND  BETAST  EQUALS' , BST, ' RADS' 
PRINT*, '  ' 

CONS=RST*COS (BST) -HST*SIN (BST) 

SBbSIN(B) 

CBsCOS (B) 

TBsTAN(B) 

CTB=1 . /TB 
XL>bCONS/SB 

C  THIS  PORTION  OF  PROGRAM  CALCULATES  A  NEW  VALUE  OF  RW 
C  NEEDED  FOR  FORCE  CONTI^ITY  BETWEEN  PHASE  II  AND 
C  TRANSITION  PHASE 

PRINT*, 'FINAL  PHASE  II  RW  EQUALS' ,RN, ' FPS' 

9M^T/ (XN*DS*CB**2 . ) -CF* (XL-XLS) / (DS*CB) 

PRINT*, ' INITIAL  VALUE  FOR  RW  FOR  TRAMS.  PHASE  EQ',RW, 'FPS' 
1700  SB»S1N(B) 

CB-COS (B) 

TBsTAN(B) 
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CTB-l./TB 

XLpCOMS/SB 

XLD«l<ir*CB 

BD>-iar*SB/xx. 

TBD>3raM 

VLXD^VCM*X0M-»-XL*CB*X0M**2 .  -XN*ar*CB*  (XL-XLS)  /XML-  (XH*OS* 

*  B]|*CB**2 . )  /XML40*SIM  (TH)  -  (Ca>L*BHO*VLT*SQRT  (VLT**2 .  +  (VCM+ 

+  «BI*XL*CB)**2.)*P*R**2.)/(2.*XML) 

K||Da«VLXD-VCN*XOM-  (XN*CB*Cr*  (XL-XLS)  ■I•XN*DS*IW*CB**2 .  -fXMC* 

*  G*SIM  (TH)  -0 . 5*C3)*RHO*  (VLT-RH)  *SQfa  {  (VLT-RN)  **2  .  +Vai**2 . ) 

*  *P*RST**2 . )  /  (XMC-t-XK*RHO*VST) 

VCMD—  (VLT-RM)  *XOM-(-  <XMC*G*COS  (TH)  -0 . 5*CD*RHO*VCM*SQRT  ( 

+  (VLT-RH) **2 . +VCM**2 . ) *P*RST**2 . ) / (XMC+XK*RHO*VST) 

XOMDo  (-VOID-  (VLT+RH)  *X0M+G*C0S  (TH)  -  (0 . 5*CDL*RHO*  (VCM+XOM* 

*  XL*CB) *SQRT (VLT**2 . + (VCN+XOM*XL*CB) **2.) *P*RP**2 . ) /XML) 

*  / (XL*CB) 

r>XM*CB*  (Cr*  (XL-XLS)  •i■DS*XLD) 

1800  xr(T-TPR)  1810,1830,1830 
1810  COMTIMDX 
1820  GO  TO  1860 
1830  PRINT*, '  ' 

PRINT*, 'TIME (SECS)','  BETA (DEG) '  ,'  FORCE (LBS)', 

+  '  VLT(FPS)','  LENGTH (FT)' 

B*57.296*B 

PRINT*, T, B, F,VLT, XL 

Bi«/57.296 

PRINT*, 'THETA (DEG) ',  '  RN(FPS)  '  ,  '  VCN(FPS)', 

*  '  OMEGA (HD/SEC)' 

TH«57.296*TH 

PRINT*, TH, RN, VCN, XOM 
PRINT*,'  ' 

1850  TPRbTPR+DTPR 
TH«TH/57.296 
I860  B-B+BD*DT 

VLTi*VLT+VLTD*DT 

TH»TH+THD*DT 

RN*iRN-»-RND*DT 

VCNWVCN+VCND*DT 

XOMbXOM-I-XOMD*OT 

ToTWT 

XF(T.6B.6.0)  GO  TO  2500 
no  TO  1700 
2500  CONTINOE 
STOP 
END 
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Appendix  B 

EXAMPLE  FUDROP  PRCXSRaM  OUTPUT 
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THB  DEFAULT  INPUTS  FOR  THIS  PROGRAM  ARB  AS  FOLLOWS 
PARACHUTE  DIAMETER  IS  28.2200FT 
CONSTRUCTED  DZA  OF  PARACHUTE  IS  28.4500FT 
MO.  OF  SUSPENSZ^  LINES  IS  28.0000 
SUSPENSION  LINE  LENGTH  IS  28.2200FT 
APERTURE  RADIUS  IS  0.820000FT 
EFFECTIVE  POROSITY  IS  l.OOOOOE-01 

SPRING  CONSTANT  OF  A  SUSPENSION  LIME  IS  20.6000LB/FT 

DAMPING  COEF.  OF  A  SUSP.  LINE  IS  1.37000LB-SBC/FT 

FIRST  PHASE  PARA.  DRAG  COEFF  IS  0.500000 

SECOND  PHASE  PARA.  DRAG  COEFF  IS  1.31500 

TOTAL  CANOPY  C  PAYLOAD  WEIGHT  IS  296.000LBS 

APP.  MASS  PROPORTIONALITY  CONSTANT  IS  2.00000 

FIRST  PRASE  CONTRACTION  FACTOR  IS  0.490000 

INITIAL  SEMI-APERTURE  ANGLE  IS  4.00000E-02RADS 

PAYLOAD  WEIGHT  IS  280.000LBS 

PAYLOAD  RADIUS  IS  0.755000FT 

PAYLOAD  DRAG  COEFF  IS  0.640000 

PARACHUTE  SNATCH  CONDITIONS  ARE  AS  FOX^WS: 

SNATCH  VELOCITY  IS  185 . OOOFT/SEC 
INITIAL  FLIGHT  PATH  ANGLE  THETA  IS  0.872700RADS 
X  COORD  OF  PAYLOAD  AT  SNATCH  IS  O.FT 
ALTITUDE  OF  PAYLOAD  AT  SNATCH  IS  lOOO.OOOFT 
OPENING  SHOCK  OF  PAYLOAD  AT  SNATCH  IS  662.000LBS 
REL.  VEL.  BETWEEN  CANOPY  AMD  LOAD  IS  O.FPS 
INITIAL  ANG.  ROTATION  RATE  IS  0. RAD/SEC 

TIME  RATE  OF  CHANGE  OF  VOLUME  IS  399.008CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  3.09262FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  80.9260FT/S 
PAYLOAD  TANCTNTIAL  ACC.  IS  -39.2486FPS/S 
ACCEL  OF  CANOPY  WRT  PLD  IS  -194 .357FPS/S 
CANOPY  NORMAL  ACC.  IS  18.3940FPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  510.826LBS 

TIME  (SECS)  VOL(COFT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

l.OOOOOE-02  7.56761  1.20761  0.825014  29.3492 

THETA(DEG)  VLT(FPS)  RW(FPS)  VCN  (FPS)  OMEGA  (RAD /SEC) 

50.0023  184.523  -3.64878  0.192856  4.67598E-04 


TIME  RATE  OF  CHAMCT  OF  VOLUME  IS  399.535CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  2.91037FT/S 

TIME  RATE  OF  CHANCT  OF  B  IS  77.5262FT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -31 .2931FPS/S 
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ACCSL  or  C&HOPY  DRT  PZD  IS  -42 . 1879rPS/S 

CAHOPy  NOBMAL  ACC.  IS  16.8182rPS/S 

RISER  rORCB  ACTING  ON  PAYLOAD  IS  441.804LBS 

TIME  (SECS)  VOL(COrT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

2.00000E-02  11.5629  1.23768  1.61925  29.3062 

THETA(DEG)  VLT(rPS)  RN(rPS)  VCN  (FPS)  OMEGA  (RAD/SEC) 

50.0028  184.171  -4.80202  0.369653  1.43896E-03 


TIME  RATE  OF  CHANGE  OF  VOLOME  IS  397.765CrS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  2.76587rT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  74.0037rT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -27 . 6413rPS/S 
ACCEL  or  CANOPY  MRT  PU>  IS  21.7000rPS/S 
CANOPY  NORMAL  ACC.  IS  15.308irPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  410.197LBS 

TIME  (SECS)  VOL(CUrT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

3.00000E-02  15.5514  1.26611  2.37862  29.2572 

TRETA(DEG)  VLT(rPS)  RW(rpS)  VCN  (FPS)  OMEGA (RAD/SEC) 

50.0040  183.877  -4.88980  0.530984  2.86631E-03 


TIME  RATE  OF  CHANGE  OF  VOLUME  IS  394.902CrS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  2.6327irT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  70.2807rT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -25 . 911irPS/S 
ACCEL  OF  CANOPY  NRT  PZJ9  IS  47.6458rPS/S 
CANOPY  NORMAL  ACC.  IS  13.7209rPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  395.317LBS 

TIME  (SECS)  VOL(CUrT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

4.10000E-02  19.9129  1.29585  3.17383  29.2051 

THETA (DEG)  VLT(rPS)  RN(rPS)  VCN  (FPS)  OMEGA (RAD/SEC) 

50.0064  183.583  -4.49456  0.691370  4.90619E-03 


TIME  RATE  OF  CHANGE  OF  VOLOME  IS  392.139CrS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  2.52727PT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  67.1337rT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -25 . 3192rPS/S 

ACCEL  or  CANOPY  NRT  PLD  IS  53.6420rPS/S 
CANOPY  NORMAL  ACC.  IS  12.3388rPS/S 
RISER  FORCE  ACTING  C»f  PAYLOAD  IS  390.331LBS 


30 


LKM6TH (FT) 


TZMB  (SBCS)  VOL(a)FT)  RADIUS (FT)  HEIGHT (FT) 

5.10000H-02  23.8495  1.32169  3.86228  29.1623 

THBTA(DHG)  VX.T(F»S)  RN(FFS)  VCN  (FPS)  OMEGA (RAD/SKC) 

50.0097  183.327  >3.98250  0.822314  7.13979E-03 


TIME  RATE  OF  CHANGE  OF  VOLOME  IS  389.489CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  2.43277FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  64.2371FT/S 
FATLOAD  TANCTNTIAL  ACC.  IS  -25 .1769FFS/S 
ACCEL  OF  CANOPY  NRT  PLD  IS  52.4640FPS/S 
CANOPY  NORMAL  ACC.  IS  11.0085FPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  389.269LBS 

TIME  (SECS)  VOL(CUFT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

6.10000B-02  27.7588  1.34653  4.52037  29.1249 

THETA(DEG)  VLT(FPS)  Rir(FPS)  VCN  (FPS)  OMEGA  (RAD/SBC) 

50.0145  183.075  -3.44765  0.939675  9.69320E-03 


TIME  RATE  OF  CHANGE  OF  VOLUME  IS  387.047CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  2.34695FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  61.5851FT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -25 .2596FPS/S 
ACCEL  OF  CANOPY  NRT  PUD  IS  48.2936FPS/S 
CANOPY  NORMAL  ACC.  IS  9.72485FPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  390.181LBS 

TIME  (SECS)  VOL(CUFT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

7.10000B-02  31.6425  1.37046  5.15061  29.0927 

THETA(DBG)  VLT(FPS)  RN(FPS)  VCN  (FPS)  OMEGA (RAD/SEC) 

50.0207  182.823  -2.94032  1.04395  1.25305E-02 


TIME  RATE  OF  CHANGE  OF  VOLUME  IS  384.837CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  2.26833FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  59.1605FT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -25 . 4541FPS/S 

ACCSL  OF  CANOPY  NRT  PLD  IS  43.0858FPS/S 
CANOPY  NORMAL  ACC.  IS  8.48i77FPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  392.082LBS 

TIME  (SECS)  VOL(CUFT)  RADIUS (FT)  HEI(fflT(FT)  LENGTH (FT) 

8.10000E-02  35.5028  1.39357  5.75537  29.0654 

THETA(DEG)  VLT(FPS)  RN(FPS)  VCN  (FPS)  OMEGA  (RAD/SEC) 

50.0287  182.570  -2.48043  1.13558  1.56191E-02 


TZMB  ItATK  or  CHAMOB  OT  VOLUME  ZS  383.043CrS/S 

TIME  KATE  OT  CHANGE  OF  RADIUS  IS  2.20285rT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  57.155€rT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -25 . 6750rPS/S 
ACCEL  or  CANOPY  NRT  PLD  IS  38.2766rPS/S 
CANOPY  NORMAL  ACC.  IS  7.40056rPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  394.208LBS 

TIME  (SECS)  VOL(CUrT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

9.00000E-02  38.9591  1.41372  £. 27988  29.0446 

THETA(DEG)  VLT(rPS)  RN(rPS)  VCN  (FPS)  OMEGA  (RAD/SEC) 

50.0374  182.340  -2.11197  1.20758  1.85896E-02 


TIME  RATE  OF  CHANGE  OF  VOLUME  ZS  381.249CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  ZS  2.13520FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  55.1055FT/S 
PAYLOAD  TANGENTIAL  ACC.  ZS  -25 . 9419FPS/S 
ACCEL  or  CANOPY  MRT  PLD  ZS  33.1880FPS/S 
CANOPY  NORMAL  ACC.  IS  6.2322irPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  396.776LBS 

TIME  (SECS)  VOL(CUrT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

l.OOOOOE-01  42.7813  1.43544  €.84186  29.0251 

THETA(DEG)  VLT(rPS)  RN(rPS)  VCN  (FPS)  OMEGA (RAD /SEC) 

50.0490  182.082  -1.75242  1.27629  2.20776E-02 


TIME  RATE  OF  CHANGE  OF  VOLUME  ZS  379.642CrS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  ZS  2.07231FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  53.223irT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -26 .2143FPS/S 
ACCEL  or  CANOPY  NRT  PLD  IS  28.5380rPS/S 
CANOPY  NORMAL  ACC.  ZS  5.09906rPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  399.410LBS 

TIME  (SECS)  VOL(CUrT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

0.110000  46.5864  1.45651  7.38432  29.0090 

THETA (DEG)  VLT(rPS)  RN(rPS)  VCN  (FPS)  OMEGA (RAD/SEC) 

50.0625  181.821  -1.44187  1.33349  2.57392E-02 


TIME  RATE  OF  CHANGE  OF  VOLUME  IS  378.201CrS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  ZS  2.01362rT/S 


TXHB  RXn  or  CHAMGS  OF  B  IS  51.4902rT/S 
PAYLOAD  TAMGKMTZAL  ACC.  IS  -26 . 4815rPS/S 
ACCBL  or  CANOPY  NRT  PLD  IS  24.3893rPS/S 
CANOPY  NOBMAL  ACC.  IS  3.99957rPS/S 
RISKR  rORCB  ACTING  ON  PAYLOAD  IS  402.019LBS 


TIMX  (SICS) 
0.120000 
TBXTA(DBG) 
50.0783 


VOL (COTT)  RADIUS (FT) 

50.3762  1.47696  7.90863 

VLT(rPS)  RN(rPS) 

181.558  -1.17557  1.37950 


HKIGHT(rT)  LKRGTa(rT) 

28.9958 

VCN  (FPS)  OMEGA (RAD/SBC) 
2 . 95522K-02 


TIME  RATE  OF  CHANGE  OF  VOLUME  IS  376.902CrS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  1.9587irT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  49.8908FT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -26 .7384rPS/S 
ACCEL  or  CANOPY  NRT  PLD  IS  20.7359PPS/S 
CANOPY  NORMAL  ACC.  IS  2 . 93252rPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  404.557LBS 


TIME  (SECS) 
0.130000 
THETA (DEG) 
50.0962 


VOL (CUFT)  RADIUS (FT) 

54.1522  1.49685  8.41623 

VLT(rPS)  RN(rPS) 

181.292  -0.948496  1.41467 


HEIGHT (FT)  LENGTH (FT) 

28.9851 

VCN  (FPS)  OMEGA (RAD/SEC) 

3.34964E-02 


TIME  RATE  OF  CHANGE  OF  VOLUME  IS  375.725CrS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  1.90717rT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  48.4103rT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -26 . 981irPS/S 

ACCEL  or  CANOPY  NRT  PID  IS  17.5555rPS/S 
CANOPY  NORMAL  ACC.  IS  1.89695rPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  406.991LBS 


TIME  (SECS) 
0.140000 
THETA (DEG) 
50.1164 


VOL (CUFT)  RADIUS (FT) 

57.9159  1.51620  8.90839 

VLT(rPS)  RN(rPS) 

181.023  -0.755820  1.43931 


HEIGHT (FT)  LENGTH (FT) 

28.9765 

VCN  (FPS)  OMEGA (RAD/SEC) 

3.75531B-02 


TIME  RATE  OF  CHANGE  OF  VOLUME  IS  374.653CrS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  1.85870rT/S 

TIME  RATE  OF  CHANCT  OF  H  IS  47.0362rT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -27 .2088rPS/S 
ACCEL  (V  CANOPY  NRT  PLD  IS  14.797irPS/S 
CANOPY  NORMAL  ACC.  IS  0 . 892078rPS/S 
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RISER  FORCE  ACTING  COI  PAYLOAD  IS 


409.312LBS 


TIME  (SECS) 
0.150000 
THETA  (DEG) 
50.1390 


VOL (COFT)  RADIUS (FT) 

61.6682  1.53505  9.38622 

VLT(FPS)  RW(FPS) 

180.752  -0.593041  1.45373 


HEIGHT (FT)  LENGTH (FT) 

28.9697 

VCN  (FPS)  OHEGA  (RAD/SEC) 

4.17051E-02 


TIME  RATE  OF  CHANGE  OF  VOLUME  IS  373.669CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  1.81300FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  45.7572FT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -27 .4205FPS/S 
ACCEL  OF  CANOPY  NRT  PLO  IS  12.4195FPS/S 
CANOPY  NORMAL  ACC.  IS  -8 .27179E-02FPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  411.510LBS 


TIME  (SECS) 
0.160000 
THETA (DEG) 
50.1640 


VOL (COPT)  RADIUS (FT) 

65.4102  1.55343  9.85075 

VLT(FPS)  RH(FPS) 

180.479  -0.456084  1.45824 


HEIGHT (FT)  LENGTH (FT) 

28.9644 

VCN  (FPS)  OMEGA (RAD/SEC) 

4.59368E-02 


TIME  RATE  OF  CHANGE  OF  VOLUME  IS  372.760CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  1.76983FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  44.5637FT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -27 . 6166FPS/S 
ACCEL  OF  CANOPY  NRT  PID  IS  10 .37396FPS/S 
CANOPY  NORMAL  ACC.  IS  -1 . 02798FPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  413.589LBS 


TIME  (SECS) 
0.170000 
THETA (DEG) 
50.1914 


VOL (COFT)  RADIUS (FT) 

69.1428  1.57137  10.30289 

VLT(FPS)  RN(FPS) 

180.204  -0.341355  1.45314 


HEIGHT (FT)  LENGTH (FT) 

28.9603 

VCN  (FPS)  OMEGA (RAD/SEC) 

5.02336E-02 


TIME  RATE  jF  CHANGE  OF  VOLUME  IS  371.914CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  1.72896FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  43.4470FT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -27 .7969FPS/S 
ACCEL  OF  CANOPY  NRT  PXD  IS  8.62326FPS/S 
CANOPY  NORMAL  ACC.  IS  -1 . 94417FPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  415.545LBS 

TIME  (SECS)  VOL(CUFT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

0.180000  72.8665  1.58888  10.7434  28.9574 
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TRXTX(DEG)  VLT(rPS)  VMiTVS)  VCM  (ITS)  0IIBGX(IUUD/8EC) 

50.2213  179.927  -0.245714  1.43871  5.45820K-02 


T2KB  RATE  OF  CSMUSX  OF  VOLUME  IS  371.122CFS/S 

TIME  RAVE  OF  CHANGE  OF  RADIUS  IS  1.69022FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  42.3995FT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -27 . 9618FPS/S 
ACCEL  OF  CANOPY  NRT  PU)  IS  7.13040FPS/S 
CANOPY  NORMAL  ACC.  IS  -2 . 83169FPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  417.381LBS 

TIME  (SECS)  VOL(CUFT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

0.190000  78.5821  1.80599  11.1731  28.9553 

THETA (DEG)  VLT(FPS)  RH(FPS)  VCN  (FPS)  OMEGA (RAD/SEC) 

50.2537  179.848  -0.188398  1.41525  5.89898E-02 


TIME  RATE  OF  CHANGE  OF  VOLXSffi  IS  370.373CFS/S 

TIME  RATE  OF  CHANGE  OF  RADIUS  IS  1.85342FT/S 

TIME  RATE  OF  CHANGE  OF  H  IS  41.4148FT/S 
PAYLOAD  TANGENTIAL  ACC.  IS  -28 .1125FPS/S 
ACCEL  OF  CANOPY  NRT  PLD  IS  5.85887FPS/S 
CANOPY  NORMAL  ACC.  IS  -3 . 89090FPS/S 
RISER  FORCE  ACTING  ON  PAYLOAD  IS  419.108LBS 

TIME  (SECS)  VOL(CUFT)  RADIUS (FT)  HEIGHT (FT)  LENGTH (FT) 

0.200000  80.2899  1.82273  11.5927  28.9539 

THETA(DEG)  VLT(FPS)  RN(FPS)  VCN  (FPS)  OMEGA (RAD /SEC) 

50.2887  179.388  -l.OlOOOE-01  1.38305  8.33854E-02 

END  OF  FIRST  FILLING  PHASE 

PARACHUTE  DYNAMICS  DURING  SECOND  FILLING  PHASE 
ITERATION  SEQUENCE  FOR  RN  BETNEEN  PHASES  I  AND  II 
FINAL  PHASE  I  BN  EQUALS  -9 .51428E-02FPS 
INITIAL  PHASE  II  RN  EQUALS  -0.129291FPS 

TINE(SECS)  BETA(DEG)  VOL(CUFT)  RADIUS (FT)  FORCE(LBS) 

0.251000  2.38848  88.9940  1.88810  488.851 

THETA(DXG)  VLT(FPS)  RN(FPS)  VCN(FPS)  OMEGA (RAD/SEC) 

50.5087  177.710  0.890433  1.08485  8.81435E-02 


TIME (SECS)  BETA (DEG)  VOL(CDFT)  RADIUS (FT)  FORCE (LBS) 

0.302000  2.50958  95.7100  1.77208  550.802 
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THBTA(DB6) 

50.7894 


VLT(FPS) 

175.723 


Rir(FPS)  VCM(FPS) 

1.55725  0.624452  0.107459 


OMEGA (RAD/SEC) 


TIME (SECS) 
0.351999 
THETA (DEG) 
51.1238 


BETA (DEG)  VOL(aJFT)  RADIUS (FT)  FORCE (LBS) 

2.66439  107.260  1.87782  629.700 

VLT(FPS)  RN(FPS)  VCN(FPS)  OMEGA (RAD/SEC) 

173.383  2.15085  6.96175E-02  0.125828 


TIME (SECS) 
0.401998 
THETA  (DBG) 
51.5060 


BETA (DBG)  VOL(CUFT) 

2.86154  122.839  2.01220 

VLT(FPS)  Rir(FPS) 

170.548  2.80460  -0.524261 


RADIUS (FT)  FORCE (LBS) 

729.786 

VCM(FPS)  C»<EGA  (RAD/SEC) 
0.140606 


TIME (SECS) 
0.451998 
THETA (DEG) 
51.9244 


BETA  (DEG)  VOL(CUFT) 

3.11067  143.852  2.18129 

VLT(FPS)  RH(FPS) 

167.084  3.55711  -1.09000 


RADIUS (FT)  FORCE (LBS) 

855.361 

VCM(FPS)  OMEGA (RAD/SEC) 

0.150909 


TIME (SECS) 
0.501997 
THETA (DEG) 
52.3655 


BETA (DBG)  VOL (CUFT) 

3.42306  172.196  2.39203 

VLT(FPS)  RW(FPS) 

162.835  4.40571  -1.56644 


RADIUS (FT)  FORCE (LBS) 

1009.81 

VCN(FPS)  (MEGA (RAD/SEC) 

0.156301 


TIME (SECS) 
0.551996 
THETA (DEG) 
52.8153 


BETA (DEG)  VOL (CUFT) 

3.81237  210.422  2.65241 

VLT(FPS)  RN(FPS) 

157.629  5.31289  -1.91017 


RADIUS (FT)  FORCE (LBS) 

1194.91 

VCN(FPS)  OMEGA (RAD/SEC) 

0.157011 


TIME  (SECS) 
0.601996 
THETA (DBG) 
53.2616 


BETA (DEG)  VOL (CUFT) 

4.29595  261.971  2.97193 

VLT(FPS)  R*r(FPS) 

151.293  6.20537  -2.10428 


RADIUS (FT)  FORCE (LBS) 

1409.94 

VCM(FPS)  OMEGA (RAD/SEC) 

0.153996 


TIME (SECS) 
0.651995 
THETA (DEG) 
53.6957 


BETA (DEG)  VOL (CUFT) 

4.89744  331.512  3.36247 

VLT(FPS)  RN(FPS) 

143.667  6.97251  -2.15939 


RADIUS (FT)  FORCE (LBS) 

1650.95 

VCN(FPS)  OMEGA (RAD/SEC) 

0.148757 
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TIME (SECS) 
0.701995 

THETA (DEG) 
54.1137 

BETA (DEG) 

5.65181 

VLT (FPS) 
134.622 

VOL(CDFT) 

425.468  3.84009 

RN(FPS) 

7.47378  -2.10603 

RADIUS (FT)  FORCE (LBS ) 

1910.76 

VCN(FPS)  OMEGA  (RAD/SEC) 

0.142974 

TIME (SECS) 
0.751994 
THETA  (DEG) 
54.5161 

BETA (DBG) 
6.61641 
VLT (FPS) 
124.076 

VOL (CUFT) 

552.986  4.42906 

RN(FPS) 

7.55701  -1.98182 

RADIUS (FT)  FORCE (LBS) 

2180 .92 

VCN(FPS)  OMEGA  (RAD/SEC) 

0.138108 

TIME  (SECS) 
0.801993 
THETA  (DBG) 

54.9071 

BETA (DBG) 
7.89872 

VLT (FPS) 

111.987 

VOL (CUFT) 

728.038  5.17214 

RN(FPS) 

7.08338  -1.81905 

RADIUS (FT)  FORCE (LBS) 

2458.30 

VCM(FPS)  OMEGA  (RAD/SBC) 

0.135155 

TIME  (SECS) 

0.851993 

THETA (DBG) 

55.2928 

BETA (DEG) 

9.75002 
VLT (FPS) 

98.2653 

VOL (CUFT) 

975.311  6.16577 

RN(FPS) 

5.93602  -1.63515 

RADIUS (FT)  FORCE (LBS) 

2770.45 

VCN(FPS)  OMEGA (RAD/SEC) 

0.134553 

TIME  (SECS) 
0.901992 

THETA (DEG) 
55.6820 

BETA (DEG) 
12.2019 
VLT (FPS) 
83.0413 

VOL (CUFT) 

1262.00  7.35082 

RN(FFS) 

0.275904  -1.51679 

RADIUS (FT)  FORCE (LBS) 

2826.11 

VCN(FPS)  OMEGA (RAD/SEC) 

0.139136 

END  OF  SECOND 

FILLING  PHASE 

BEGINNING  OF 

TRANSITION 

PHASE 

STEADY  STATE  CANOPY  VOLUME  EQUALS  1292.47CU  FT 

STEADY  STATE  CANOPY  RADIUS  EQUALS  7.48724rT 

HST  EQUALS  0.829277rT.  AND  BBTAST  EQUALS  0.218339RADS 

FINAL  PHASE  II  RM  EQUALS  -1.55267FPS 

INITIAL  VALUE  FOR  RN  FOR  TRANS.  PHASE  EQ  3.33678FPS 


TIMB<SECS) 
1.01099 
THETA  (DBG) 
58.6914 


BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

13.0902  976.350  55.4796  31.4806 

Rir(FPS)  VCN(FPS)  OMBGA(RD/SEC) 

-23.5077  -2.10969  0.190587 
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TIME (SECS)  BETA (DEG)  FORCE (LBS)  VLT(rPS)  LENGTH (FT) 

1.11000  14.0052  122.299  48.9713  29.4811 

THETA(DEG)  RN(FPS)  VCN(FPS)  OMEGA(RD/SEC) 

57.9422  -15.8475  -2.82439  0.250439 


TIME (SECS)  BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

1.21000  14.4675  42.6589  47.8726  28.5388 

THETA (DEG)  RN(FPS)  VCN(FFS)  CBIEGA (RD/SEC) 

59.5186  -3.76408  -3.52420  0.296657 


TIME (SECS)  BETA(DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

1.31001  14.4890  217.216  46.0274  28.4974 

THETA (DBG)  RN(FPS)  VCN(FPS)  OMEGA (RD/SEC) 

61.3012  1.73288  -4.12543  0.323326 


TIME (SECS)  BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

1.41001  14.3900  321.982  42.8341  28.6893 

THETA  (DEG)  RN(FPS)  VCN(FPS)  CBffiGA  (RD/SEC) 

63.2045  1.66064  -4.64702  0.340446 


TIME (SECS)  BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

1.51001  14.3461  314.612  39.6048  28.7752 

THETA(DBG)  RN(FPS)  VCN(FPS)  OHEGA(RD/SBC) 

65.1974  0.121788  -5.12113  0.355167 


TIME (SECS)  BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

1.61002  14.3632  268.510  37.1501  28.7417 

THETA(DEG)  RN(FPS)  VCN(FPS)  ^CEGA(RD/SEC) 

67.2703  -0.638844  -5.55814  0.368046 


TIME (SECS)  BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

1.71002  14.3948  237.452  35.4439  28.6800 

THETA(DBG)  RN(FPS)  VCN(FPS)  QMEGA(RD/SEC) 

69.4071  -0.543065  -5.95372  0.377165 


TIME (SECS)  BETA(DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

1.81003  14.4137  228.468  34.1564  28.6432 

THETA(DEG)  Rir(FPS)  VCN(FPS)  OMEGA(RD/SEC) 

71.5823  -0.220812  -6.30131  0.381293 
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Txia(sscs) 

1.91003 

TBBTA(DK6) 

73.7668 


TZMB (SICS) 
2.01004 
TBI7A(PIG) 
75.9335 


TIME (SICS) 
2.11003 
THITA(DBG) 
78.0580 


BKIAK 


TIMB(SECS) 

2.91097 

THBTA(DEG) 

91.1987 


TIMB(SECS) 
3.01096 
THETA  (DIG) 
92.1497 


TIME (SECS) 
3.11096 
THETA  (DBG) 
92.9256 


TZMB (SECS) 
3.21095 
THETA (DEG) 
93.5271 


BETA  (DEG)  rOBCE(ZaS)  VLT(rES)  LENGTH (FT) 

14.4195  228.304  33.0598  28.6317 

EM  (EPS)  VCH(rPS)  OMEGA  (RD/SEC) 

4. 755311-02  -6.59731  0.380441 


BETA (DEG)  rOBCE(LBS)  VLT(rPS)  LENGTH (FT) 

14.4212  226.931  32.0906  28.6286 

RM(FPS)  VCN(FPS)  OMBGA(BD/SEC) 

-3.62383E-02  -6.84028  0.375116 


BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4239  222.614  31.2553  28.6233 

RM(FPS)  VCM(FFS)  OMEGA (RD/SBC) 

-7.37384E-02  -7.02922  0.365730 


IN  PRZNyODT  HERE  TO  SAVE  SPACE 


BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4518  193.173  27.8949  28.5692 

B»(FPS)  VCN(FPS)  OMBGA(RD/SEC) 

-5.22438E-02  -6.55442  0.181046 


BETA  (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4543  190.705  27.7164  28.5644 

RN(FPS)  VCN(FPS)  omega (RD/SEC) 

4.64431E-02  -6.26577  0.150566 


BETA (DBG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4564  188.582  27.5726  28.5602 

R»(FPS)  VCN(FPS)  OMEGA (RD/SBC) 

-4.02695E-02  -5.93584  0.119985 


BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4583  186.808  27.4599  28.5566 

RN(FPS)  VCN(FPS)  OMEGA  (RD/SBC) 

-3.38517E-02  -5.56892  8. 975521-02 


39 


TZME(SBCS)  BBTA(DBG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

3.31094  14.4598  185.387  27.3752  28.5537 

TRETA(DEG)  RN(FFS)  VCN(FPS)  OMBGA(RD/SEC) 

93.9573  -2.72760E-02  -5.16971  6.03079E-02 


TIME (SECS)  BETA(DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

3.41094  14.4610  184.320  27.3150  28.5513 

THETA  (DEG)  Rir(FPS)  VCN(FPS)  OMEGA  (RD/SEC) 

94.2220  -2.07229E-02  -4.74327  3.20465E-02 


TIME  (SECS)  BETA(DEG)  FORCE  (LBS)  VLT(FPS)  XJENGTH(FT) 

3.51093  14.4619  183.587  27.2762  28.5496 

THETA  (DEG)  BK(FPS)  VCN(FPS)  C3MEGA(RD/SEC) 

94.3291  -1.44942E-02  -4.29498  5.34339E-03 


TIME (SECS)  BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

3.61092  14.4625  183.153  27.2552  28.5485 

THETA  (DEG)  RN(FPS)  VCN(FPS)  OMEGA  (RD/ SBC) 

94.2884  -8.89619E-03  -3.83049  -1.94656E-02 


TIME (SECS)  BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

3.71091  14.4628  182.989  27.2486  28.5479 

THETA (DEG)  RN(FPS)  VCN(FFS)  OMEGA (RD/ SEC) 

94,1116  -4,01141E-03  -3.35568  -4.20845E-02 


TIME (SECS)  BETA(DEG)  FORCE (LBS)  VLT(FFS)  LENGTH (FT) 

3.81091  14.4629  183.026  27.2532  28.5477 

THETA  (DEG)  RN(FPS)  VCN(FPS)  OMEGA  (RD/SBC) 

93.8120  -8.20918E-05  -2.87656  -6.22611E-02 


TIME(SECS)  BETA(DE6)  FORCE(LBS)  VLT(FPS)  LEMGTH(FT) 

3.91090  14.4628  183,219  27.2657  28.5479 

THETA  (DEG)  RN(FPS)  VCN(FPS)  CMEGA  (RD/SEC) 

93.4043  2.92303E-03  -2.39922  -7.97895E-02 


TIME (SECS)  BETA (DEG)  FORCE (IAS)  VLT(FFS)  UEMGTR(FT) 

4.01089  14.4627  183.474  27.2831  28.5482 
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THETA  (DBG) 
92.9040 


TIME (SECS) 
4.11088 
THETA (DBG) 
92.3277 


TZMB (SECS) 
4.21088 
THETA  (DEG) 
91.4919 


TIME (SECS) 
4.31087 
THETA (DBG) 
91.0139 


TIME (SECS) 
4.41086 
THETA (DEG) 
90.3105 


TIME (SECS) 
4.51086 
THETA (DBG) 
89.5982 


TZMB (SECS) 
4.61085 
THETA  (DBG) 
88.8926 


TIME (SECS) 
4.71084 
THETA (DBG) 
88.2088 


EN(rPS)  VCM(rPS)  0MBGA(HD/SEC) 

4.92302E-03  >1.92973  -9.45134E-02 


BETA(DBG)  rOHCE(LBS)  VLT(rPS)  LENGTH (ET) 

14.4623  183.860  27.3025  28.5488 

HN(rPS)  VCN(rPS)  0MBGA(BD/SEC) 

5.84737B-03  -1.47400  -0.106327 


BETA(DBG)  rOHCE(LBS)  VLT(rPS)  LENGTH (ET) 

14.4620  184.205  27.3214  28.5495 

EM  (EPS)  VCM(EPS)  OMEGA  (RD/SBC) 

5.22327E-03  -1.03772  -0.115176 


BETA (DBG)  EOHCB(LBS)  VLT(FPS)  LENGTH (FT) 

14.4618  184.417  27.3387  28.5499 

RN(FPS)  VCN(EPS)  OMEGA  (RD/SBC) 

5.14782E-03  -0.626208  -0.121060 


BETA (DEG)  FORCE (LBS)  VLT(EPS)  LENGTH (FT) 

14.4616  184.654  27.3526  28.5503 

Rll(FPS)  VCN(EPS)  OMEGA  (RD/SEC) 

4.81418B-03  -0.244295  -0.124025 


BETA  (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4614  184.813  27.3620  28.5506 

RM(FPS)  VCN(rPS)  OMEGA (RD/SEC) 

4.09617E-03  1.03768E-01  -0.124171 


BETA  (DBG)  FORCE  (LBS)  VLT(FPS)  UMGTB(FT) 

14.4612  184.918  27.3660  28.5510 

Rir(FPS)  VCN(FPS)  OMEGA  (RD/SBC) 

1.88966B-03  0.414419  -0.121641 


BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4612  184.835  27.3649  28.5510 

Rir(FPS)  VCN(FPS)  OMEGA  (RD/SEC) 

-7.62786B-04  0.684843  -0.116623 
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TZMB(SXCS)  BXTA(DEG)  rORCB(LBS)  VLT(FPS)  LENGTH (FT) 

4.81083  14.4613  184.645  27.3594  28.5508 

IW(FPS)  VCN(rFS)  OHEGX(HD/SBC) 

87.5602  -2.82502E-03  0.913006  -0.109339 


TIME (SECS)  BETA(DEG)  FORCE (LBS)  VLT(FPS)  LENGTH(FT) 

4.91083  14.4615  184.424  27.3504  28.5504 

THETA(DEG)  RW(FPS)  VCN(rPS)  OMEGA(RD/SEC) 

86.9592  -3.85151E-03  1.09768  -1.00042E-01 


TIME (SECS)  BETA(DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

5.01082  14.4617  184.164  27.3382  28.5501 

THETA  (DEG)  RM(FPS)  VCN(FPS)  C»ffiGA(RD/SEC) 

86.4165  -5.31595E-03  1.23846  -8.90130E-02 


TIME (SECS)  BETA(DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

5.11081  14.4620  183.792  27.3246  28.5494 

THETA (DEG)  RN(FPS)  VCN(FPS)  OMEGA (RD/SEC) 

85.9412  -5.45431E-03  1.33571  -7.65486E-02 


TIME (SECS)  BBTA(DBG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

5.21080  14.4623  183.539  27.3102  28.5489 

THETA(DEG)  RN(FPS)  VCN(FPS)  OMEGA(RD/SEC) 

85.5407  -5.11840E-03  1.39052  -6.29591E-02 


TIME (SECS)  BETA(DBG}  FORCE (LBS)  VLT(PPS)  LENGTH (FT) 

5.31080  14.4626  183.222  27.2959  28.5484 

THETA (DEG)  RN(FPS)  VCN(FPS)  OMEGA (RD/SEC) 

85.2205  -5.00715E-03  1.40468  -4.85611E-02 


TIME (SECS)  BETA(DBG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

5.41079  14.4628  182.973  27.2826  28.5479 

THETA  (DEG)  RN(FPS)  VCN(FPS)  OffiGA  (RD/SEC) 

84.9844  -5.01031E-03  1.38058  -3.36704E-02 


TIME (SECS)  BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

5.51078  14.4630  182.775  27.2709  28.5476 

THETA (DEG)  RN(FPS)  VCN(FPS)  OMEGA (RD/ SEC) 

84.8342  -5.08950E-03  1.32112  -1.85966E-02 
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TIMB(SBCS) 

5.61078 

TBBTA(DIG} 

84.7702 


TIMBCSKCS) 

d  5.71077 

THETA  (DBG) 
84.7909 


TIME (SECS) 
5.81076 
THETA (DEG) 
84.8933 


TIME (SECS) 
5.91075 
THETA (DBG) 
85.0729 


BBTA(DEG)  FORCE (LBS)  VLT(rPS)  LENGTH (FT) 

14.4632  182.611  27.2619  28.5472 

Rir(FPS)  VCN(FFS)  OMEGA  (HD/SEC) 

-4.51721E-03  1.22970  -3.63787E-03 


BETA (DBG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4633  182.489  27.2561  28.5469 

RW(FPS)  VCN(FPS)  OMEGA(RD/SBC) 

-2.80738B-03  1.11005  1.09237E-02 


BETA (DEG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4634  182.502  27.2535  28.5468 

HW(FPS)  VCN(FFS)  OMEGA  (RD/SEC) 

-3.17685E-04  0.966224  2.48261E-02 


BETA (DBG)  FORCE (LBS)  VLT(FPS)  LENGTH (FT) 

14.4634  182.561  27.2533  28.5468 


R]r(FPS) 

1.33544B-03 


VCN(FPS) 

0.802478 


OMEGA  (RD/SBC) 
3.78316E-02 
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